NA' L p 
MG: 



L 1 -Y 



/.HOOL 



A NUMERICAL METHOD OF 1000-MB PROGNOSIS 

THROUGH 

PREDICTION OF THICKNESS FIELD 






Ho-Chih Chin 



Chuang-Tao Kuan 
Ha Song Pak 



A NUMERICAL METHOD OF 1000-MB PROGNOSIS 

THROUGH 

PREDICTION OF THICKNESS FIELD 



by 



Ho-Chih Chin Commander 

Republic of China Navy 

Chuang-Tao Kuan Lieutenant Commander 

Republic of China Navy 

Ha Song Pak Lieutenant (Junior Grade) 

Republic of Korea Navy 



Submitted in partial fulfillment of 
the requirements for the degree of 

MASTER OF SCIENCE 
IN 

METEOROLOGY 



United States Naval Postgraduate School 
Monterey, California 



19 6 1 



VslV'S r-V 
\<\(*\ 



L 



• -OL 



A NUMERICAL METHOD OF 1000-MB PROGNOSIS 

THROUGH 

PREDICTION OF THICKNESS FIELD 
by 

Ho-Chih Chin 
Chuang-Tao Kuan 
Ha Song Pak 

This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 

METEOROLOGY 
from the 

United States Naval Postgraduate School 



ACKNOWLEDGEMENTS 



The writers are indebted to Professor George J. 
Haltiner of the Department of Meteorology and Oceanography 
at United States Naval Postgraduate School for suggesting 
the topic, for outlining the way in which the study might 
be pursued , and for the assistance he rendered them. 

They are also indebted to Lieutenant Commander Edward 
M. Carlstead of the Fleet Numerical Weather Facility, 
United States Navy, for the assistance he gave them during 
the course of investigation. 



ii 



ABSTRACT 



A two-level numerical prediction model incorporating 
terrain and non-adiabatic warming effects is tested. An 
attempt is made to determine the empirical coefficient (s) 
for each of three terms composing a prognostic equation 
designed to predict the 1000 to 500-mb thickness through 
the use of a digital computer. 

The model is applied to three cases in April for a 
large part of the Northern Hemisphere. The predicted 
positions of most of the pressure systems were reasonably 
accurate; however the model tends to over-develop pressure 
centers. The effects of terrain and non-adiabatic warming 
were in accord with theoretical expectations; however, the 
excessive development dominated the numerical verification 
so that no significant improvement resulted from the in- 
clusion of these refinements. 
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1. INTRODUCTION 



The availability of high-speed large-capacity elec- 
tronic digital computers has long since made numerical 
weather prediction an operational reality. Notable success 
has been achieved at the 500-mb level. However, surface 
prognoses obtained by numerical integration techniques on 
the average have not equalled the accuracy obtained by 
subjective methods. Even with some relatively complex 
models, the degree of success has not corresponded to the 
degree of sophistication of the models. To improve the 
quality of surface prognosis, two alternatives seem to be 
feasible. One approach is to devise even more elaborate 
mathematical models in order to take into account as many 
pertinent weather parameters as possible. Perhaps new 
parameters which are not just a recombination of the 
familiar physical variables — pressure, temperature and 
humidity, will have to be introduced and measured to a high 
degree of accuracy on a synoptic time and space scale. 
Another alternative is to try quasi-empirical prediction 
models with a reasonable dynamical basis. In this paper 
a simple quasi-empirical 1000-mb prognostic model is 
derived and testing results are presented. 

Since 500-mb prognoses from the barotropic model or 
its variations have been quite successful, a next logical 
step would be to devise a method to forecast the thickness 
field; and then through the predicted thickness arrive at 
the prognosis of some other level. In this experiment, a 
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technique for predicting the thickness field between 
1000 mb and 500 mb is combined with results of the baro- 
tropic model currently in use by the Fleet Numerical 
Weather Facility. The predicted 500-mb stream field at 
every hour is applied in turn to generate a thickness 
between levels through a thickness-tendency equation. 

After sufficient iterations, a 24-hour 1000-mb prognos o 
contour map is obtained. 

From the hydrostatic viewpoint, pressure is simply the 
weight of the air column above a particular level, and 
local changes in pressure may be thought of as the inte- 
grated density change in the air column. The density change 
can be described in terms of the temperature field which, 
in turn, can be expressed as a thickness change. The main 
feature of this model, the thickness-tendency equation, 
takes into account the temperature field, horizontal 
motions, and the vertical motions induced by terrain and 
surface stress. 

Actually there is no inherent restriction on the choice 
of the 500-mb prediction equation. The thickness-tendency 
equation could be combined with any 500-mb model, though 
perhaps requiring some minor modifications in the pro- 
gramming. Obviously, if another 500-mb model is employed 
instead of the particular barotropic model used here, the 
results of testing may be quite different, and the empir- 
ical coefficients may not have the same values as obtained 
in this experiment. 
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2. Prognostic Model 



The prognostic model is composed of two prediction 
equations. The first is the 500-mb stream-barotropic vor- 
ticity equation, the other is a thickness-tendency equation. 

With the twisting term and vertical advection of vor- 
ticity omitted, the vorticity equation in (x,y,p,t) coor- 
dinates becomes 

+v*v(M)= -Of+f)vv , (1) 

Due to the fact that at intermediate levels in the tropo- 
sphere between 500 and 600 mb, the horizontal velocity 
divergence is normally relatively small and frequently may 
be neglected, the vorticity equation may be written in the 
form 

i £- +V'P(f+i)= 0 . (2) 

Mere geostrophic approximation of (2) with 

V 3 - 4= a 

yields the barotropic quasi-geostrophic model 

V‘(£) = -J(*. fWf) m 

where J is the Jacobian operator. Equation (3) has had 
considerable success in predicting the height field for 
500 mb; but due to the velocity divergence inherent in the 
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geostrophic wind, an error commonly referred to as "spuri- 



ous anticyclogenesis" is introduced. To remove this source 
of error, a non-divergent wind is used as follows: 

%= -vt»t , i= V*Y . 

With these substitutions, equation (2) becomes: 

VdD+TC*, + f ) = o , <4, 

A characteristic of the barotropic model as repre- 
sented by (4) is that very long waves, say number 1, 2, and 
3 with respect to a latitude circle, tend to retrogress at 
very high speeds, a phenomenon not observed in nature. To 
control this error, an additional term has been introduced. 
The result is the Helmholtz stream-barotropic equation 




The first law of thermodynamics and the Poisson’s 
equation can be combined in the form 




Also from Poisson's equation and the equation of state, it 
follows that 

r\ , n i -f^ v / Cr -i < 7 > 

0 = d R, T lc0 ° 
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which may be used in (6) to give 



£l d- a-\/.\~7A -+ —A— £&— o = Q cl (8) 

at + VVa+ e <*f c+ t . 

Finally replacing (j^ in the first two terms on the left 
by means of the hydrostatic equation, and on the right 
through the equation of state yields 

&(%) = -V'Vdf)-^-^- ». 

Here (J — -r— r is the stability factor and is 

u 6 by df 

considered to be constant in space and time. 

The "vertical velocity" may be considered to consist 
of three contributions, the large-scale vertical motion, 
the terrain-induced vertical motion and f rictionally-in- 
duced vertical motion. Here these contributions will be 
represented by assuming 60 to be of the form 



CO — : co 5 



1- (^fl -?f.*^ J VT*vHt + w l T 



(10) 



where 



,1 



1= 



t- 






X < V ^ X 



l 



y 
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The first term on the right side of (10) is the para- 



bolic distribution frequently assumed to be typical of 
large-scale pressure systems [1] . 

The second term of CO represents the terrain-induced 
vertical motion. Here the assumption is made that the sur- 
face vertical velocity Wq is proportional to the horizontal 
wind and the gradient of terrain height at the surface; 
i.e. , Wq = k^ • H^. Since Wq, it follows that 



city is further assumed to decrease linearly with pressure; 
thus we arrive at the form in the second term. 

The third term represents vertical motion induced by 
surface friction. Here U, is the "vertical velocity" at 



top of the friction layer due to surface stress. It is 
assumed to have a simple distribution with respect to 
height which can be approximated linearly as shown in fig. 1. 




The terrain-generated vertical velo- 




L0— 0 at f s 






0 at % 



Figure 1. Profile of vertical velocity due to 
surface stress. 
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Cressman [2] deduced the following form for Ul in terms 
of geostrophic surface wind and a drag coefficient; 




1% 



(G V a /u*4 7*) 



11 

f 



F 



(id 



where 

F = 




( 12 ) 



Here u and v are the x- and y-components of the geostrophic 
wind at 1000 mb, f is the coriolis parameter and is the 
drag coefficient. The latter is mainly a function of the 
roughness of the underlying surface and stability in the 
micrometeorological layer, and has been determined for 
average conditions by Cressman for large-scale forecasting 
over the Northern Hemisphere. 

Upon integration between 'p and j ^ equation (9) 



becomes 



ah 

at 



-Vs’Vh Ht 



rjj F 

d| 



G + fT,-T t )g., 



+ aWl-A-Q 



(13) 



where Q represents the amount of heat added to the layer 
per unit time. 

If in equation (13) we 

1. assume (T&%= f'V- V h and combine it with 
the first term and affix a coefficient K-^; 
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2 . combine coefficients in the second term and 
designate the combination as 

3. assume Q is proportional to V » V T , or 

i ° s 

Q = K 3 V «^T S , where T g is the sea-surface 
temperature; and 

4. assume a reasonable pressure level as top 
of friction layer in the frictional term 
and designate the combined coefficient as 




we arrive at : 

at" = Mt - KM’V H t + W VV - K*F . 1141 



This is the thickness-tendency equation, with the empirical 
coefficients to be determined through testing. The thick- 
ness equation together with the 500-mb stream-barotropic 
equation constitute the prognostic model. 

The local rate of thickness change is thus seen to be 
due to the combined effect of thickness advection between 
layers, terrain effects, non-adiabatic heating and a 
frictional effect. Since cold advection is usually 
associated with adiabatic warming due to subsidence aloft 
and by heating from below in the surface layers, while 
warm advection is usually counteracted by cooling due to 
large-scale upward motion and in the surface layers by 

cooling from below, the coefficient might be expected 
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to have different values for cold and warm advection. This 



was borne out by the results of this numerical experiment. 

The terrain term represents vertical motion induced 
by uneven terrain, and obviously it vanishes over oceans 
and level land. 

The non-adiabatic term in the form re P resents 

vertical diffusion of heat between air and surface water, 
For an air mass moving over warmer water, the surface 
temperature quickly adjusts to that of the surface water 
and rapid vertical diffusion of heat takes place. The rate 
of diffusion will normally increase with increasing hori- 
zontal wind which enhances the turbulence, and upward dif- 
fusion of heat will certainly continue as long as the water 
temperature increases downwind. For warm air moving over 
colder water the stability created by surface cooling in 
the lowest layer will tend to inhibit diffusive heat ex- 
change. This suggests the empirical coefficient may 
also be assigned two values depending on the sign of 

W V 

The effects of surface stress about closed pressure 
systems may be summarized as follows: 

Low: upward motion at top of friction layer, AO ^ 

negative; contribution to thickness tendency, 
negative . 

High: downward motion at top of friction layer 

positive; contribution to thickness tendency, 
positive . 
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Hesselberg and Friedmann found that the friction effect 
is at least two orders of magnitude smaller than some of the 
other terms in the vorticity equation. The assumption is 
made that it will also be quite small as compared to other 
terms in the thickness-tendency equation. Hence in the 
present experiment it was neglected. Thus equation (14) 
becomes 



ah 

at 



-KM’Vh-K.Vt'VHt + K. ? Y • V 1j 



(15) 



where represents the wind at the terrain height. Next 
the winds and in equation (15) are approximated by 
geostrophic values resulting in 



lL_ _ X 

dt~j 



KiJ C^5. M-KjTC^Ts) K 2 \£"VH t < 16 > 



The wind at terrain height is computed by linear inter- 
polation. On the assumption that wind hodograph is approx- 
imately a straight line between 1000 and 500 mb, may be 
expressed in the form 







( 17 ) 
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Thus the wind at terrain height is equal to the wind at 
500-mb level minus some fraction of the thermal wind, and 
the fraction is the ratio of the surface/500-mb thickness 
to the 500/1000-mb thickness. With the geostrophic approx- 
imation for the wind and the letter A to designate 
it follows that 







(18) 



Hence the terrain term in equation (16) may be computed 
using two Jacobians. However, because of limited computer 
storage capacity, instead of the two Jacobian operations, the 
following transformation is made 




$ 



In finite differance form, the last expression becomes 



V t 5 VH t =l^|[-(^-AL) + (z„-AKs)](Mt e -Ht N ) 

4 -' V 



+ 



(£y e Alt e ) (%v"AhJ (H-tN+ks) 



( 20 ) 
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where the subscripts N, S, E and W are notations for the 
location of grid points. The quantity A has three var- 
iables, but may be simplified by taking standard atmos- 
pheric values for and h without introducing too large 
an error. With the letters Z fcN , Z t g, Z^ and Z^ w to denote 
the quantities within the brackets in the final expression 
for above, equation (15) has the form 



3 h 

at 




vT( 



vM-kJ dxi 




© 
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3. Procedure 



Three days' data in April 1955 were tested. Progno- 
ses were made for April 2nd, 3rd and 4th and were labeled 
as Cases 1, 2 and 3, respectively A CDC-1604 digital 
computer was employed for the computations over a 1977 
point octagonal grid for the Northern Hemisphere poleward 
of ION. At the end of each iteration, a predicted thick- 
ness is obtained at each grid point. This is subtracted 
from the prognostic 500-mb stream function to find the 
1000-mb height. The final output is a predicted 1000-mb 
height field in a contour map form. The forward-difference 



and a centered-dif ference method was used for succeeding 
hours. The forward and centered difference methods may be 
expressed respectively as 



A schematic block diagram of the thickness prediction 
program is given in fig. 2. 

Using the computer, a "pillow" and a RMSE for each 



method was used for calculation of 



hs 

9t: 



for the first hour 




and 
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Figure 2. Schematic block diagram showing procedure, 
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prognostic map were computed by comparison with the actual 
map for verification purposes. The formulae used in com- 
puting the pillow and RMSE are 



pillow 




RMSE 




where A and B stand for predicted and verifying values, 
respectively, at the same grid point,* and x stands for the 
total number of points differenced. Qualitatively, a pillow 
is the algebraic mean difference over the whole grid between 
the two fields A and B,* while the RMSE is a numerical meas- 
ure of the prognosis, and it reflects errors in position 
and in intensity of pressure systems. 

Due to fictitious reflections on the boundary points 
and inaccuracies in the geostrophic approximation at low 
latitudes, it was decided to limit the verification to the 
areas north of the 20th latitude circle. 

Since the major contribution to thickness change is 
due to the advection term, it was investigated first. 

Values of determined by previous investigators [3] were 
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used initially, then the value was varied to arrive at the 
best verification. Furthermore, since warm and cold ad- 
vection did not necessarily give the same contribution 
to the predicted change in thickness, combinations of two 
different values were tried during the investigation of 
this term. 

With the coefficients for the advection term tenta- 
tively determined, the investigation was carried on to 
find suitable coefficients for the terrain term. Only one 
coefficient of best fit was sought for the terrain term. 
During the investigation of this term, a smoothed terrain- 
height field is required in order to obtain the signifi- 
cant scale of vertical motion. Since the grid distance is 
about 380 km, any smaller- scale variation in the terrain 
field would not be included. For this reason, the smoothed 
terrain-height field data of J. Smagorinsky was used. 

Finally, an investigation was conducted to find suit- 
able coefficients for the non-adiabatic term. Two coef- 
ficients were sought for this term, one for heating and 
one for cooling. Values of determined by previous inves- 
tigators [ 4 J were also used initially in this case. The 
monthly-mean Sea-temperature field of April was used during 
the investigation of the non-adiabatic term. This field 
was obtained by graphical addition of the monthly-mean 
air-temperature charts of April and the 3-month-mean air- 
sea temperature difference charts [5] . The resulting 
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temperature field is shown in fig. 10. The temperature 
field on the land area was included to facilitate the 
computations; however, these values did not enter into the 
final results. 
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4. Results and Discussion: 

Figures 3 to 8 indicate the initial maps and the 
prognostic maps using different combinations of the terms 
forming the thickness-tendency equation. The isolines 
are drawn at intervals of 200 feet, and labeled in hun- 
dreds of feet. Case 3 has been chosen for illustration 
purposes because it gave the best results among the three 
cases tested, both in the prediction of the positions of 
the pressure centers and in numerical verifications. 

By comparing figs. 4 and 5, it can be seen that the 
positions of the pressure systems are in close agreement 
with four exceptions, namely: the high over the east coast 

of Asia; the high off the west coast of the United States; 
the high over the Caspian Sea; and the high lying to the 
west of Gibraltar. Possible reasons for errors in their 
locations will be discussed at the end of this section. 

It was also of interest to notice that the positions of 
systems in the Atlantic Ocean were predicted more accurately 
than those in the Pacific Ocean. 

Before discussing the coefficients found thus far in 
this investigation, it is necessary to mention two phenom- 
ena related to development. These phenomena are fictitious 
anticyclogenesis and cyclogenesis which predominated in 
all three cases tested. While no dynamical explanation 
could be determined from this investigation, it was apparent 
that the effect was due mainly to the advection term. 
Comparison of figs. 4 and 6 shows that most of the highs 
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were too high by an average amount of 200 feet, while most 
of the lows were too low by an average amount of 100 feet 
when the advection term was used alone to make the forecast. 
Since the major contribution to thickness change is due to 
the advection term, these phenomena invariably concealed 
the effects of a refinement like the non-adiabatic heat 
exchange in the numerical verification. 

Aside from the effect of over-prediction by the advec- 
tion term, it was found that the combination of coeffi- 
cients 0.4 for warm advection and 0.5 for cold advection 
gave the best verification when this term is used alone to 
make the forecast. Other sets of coefficients, such as 
0.8/0. 9, 0.7/0. 8, 0.5/0. 5, 0. 3/0.4, to mention only a few, 
gave larger RMSE values. Therefore, it was decided tenta- 
tively to use 0.4/0. 5 as the coefficients of best fit for 
the advection term. The results of numerical verification 
for prognoses using the advection term only with coeffi- 
cients of 0.4 for warm advection and 0.5 for cold advection 
are shown in Table 1. 
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Case 

Number 


Pillow 

(ft) 


RMSE 

(ft) 


1 

(0000Z April 2, 
1955 ) 


0 


+ 189 


2 

(0000Z April 3, 
1955 ) 


+ i4 

j 


4- 171 

1 

1 


3 

(0000Z April 4, 
1955 ) 


I 


168 

1 



Table 1. Results of numerical verification for prognoses 
using the thickness advection term only with 
coefficients of 0.4 for warm advection and 0.5 
for cold advection. 
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For the investigation of the terrain term, verifica- 
tion by pillow and RMSE in Table 2 showed no improvement 
in RMSE when the terrain term was included in this model. 

In this connection, it should be noted that the empirically- 
determined best-fit coefficient of the terrain term is 0.01, 
while the maximum terrain-height gradient is about 12 times 
larger than the thickness gradient. Hence the terrain term 
contribution can be at most about 25% of the advection term. 
Also, less than half of the total area covered by the grid 
point is affected by the terrain term, while the pillow 
and RMSE were computed over the entire grid. Therefore, 
the actual contribution of the terrain term should be 
larger than that shown by RMSE. 

In order to investigate the effects of the terrain- 
induced vertical motion in more detail a differential 
analysis was made and presented in fig. 9. Here the dif- 
ference between the forecast with thickness advection only 
(fig. 6) and that with both the terrain term and thickness 
advection (fig. 7) was computed. The isolines are drawn 
at intervals of 50 feet, and the figures labeled in 
tens of feet. Of the 20 closed isolines observed in fig. 9, 
14 indicate a decreased error; 5 of them an increased er- 
ror; and one, neutral. Thus, there are some areas where 
the terrain term does not give any improvement and, in 
fact, made the results worse. For example, over Greenland, 
the forecast made with the thickness advection alone over- 
intensifies the low, and with the inclusion of the terrain 
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Case 

Number 


Pillow 

(ft) 


RMSE 

(ft) 


1 

(0000Z April 2, 
1955 ) 


- 15 


+ 191 


2 

(0000Z April 3, 
1955 ) 


+ 1 


+ 172 


i 

j 

3 

(0000Z April 4, 
1955 ) 


- 7 


J 

+ 168 

i 



Table 2. Results of numerical verification for prog- 
noses using advection and terrain terms with 
coefficients of 0.4 for warm advection, 0.5 
for cold advection and 0.01 for the terrain 
term. 
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term, further deepening is observed on the leeward side of 
the mountain, giving an even poorer result. 

Fourteen closed isolines of the 20 on the differential 
analysis showed a better forecast, with a two- to eight-mb 
improvement in the surface pressure field , when the 
terrain-induced vertical motion over land areas was in- 
cluded. These results indicate some justification for the 
inclusion of the terrain term in the model. 

Since the terrain-induced vertical motion is a func- 
tion of the wind component parallel to the terrain-height 
gradient as well as the terrain-height gradient itself, 
large values of terrain contribution may be expected in 
areas where above-mentioned values are large. In this 
respect, the test data available for this investigation 
were not favorable for a good evaluation of the terrain 
term. The synoptic maps show that the areas where the 
terrain gradient is large are mainly associated with 
small wind components parallel to the terrain-height 
gradient . 

With only three days' data tested, no final conclu- 
sion can be reached regarding the terrain-induced verti- 
cal motion in this model at this time. Effects of the ter- 
rain-induced vertical motion were produced in areas where 
they were expected from dynamic considerations; however, the 
values are too small in general. This seems to suggest 
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that further experiments may indicate a larger coeffi- 
cient for the terrain term with perhaps some other addi- 
tive empirical factors. 

The results of numerical verification for prognoses 
using advection and non-adiabatic terms are shown in Table 
3. The eight blanks in Table 3 arose in cases where the 
forecast map was very much distorted, in contrast to the 
regular appearance of the remaining forecast maps. This 
erratic behavior was apparently connected with a malfunc- 
tion of the magnetic tape units. 

The coefficients of best fit for the non-adiabatic 
term, as shown in Table 4, were 10 for heating and 6 for 
cooling. They were chosen only because they provided the 
smallest pillow. The RMSE obtained by using this set of 
coefficients was -j- 168 feet for case 3, which equaled 
exactly that obtained by using advection term alone for 
prognosis. As a matter of fact, Table 3 shows that the 
RMSE values increased when the coefficients were increased 
algebraically. Since it was mentioned earlier that the 
RMSE is a numerical measure of the accuracy of the forecast, 
it appears that the addition of the non-adiabatic term at 
best made no improvement in the forecast. Two arguments 
need to be made clear at this point. Firstly, fictitious 
anticyclogenesis and cyclogenesis by the advection term 
alone already have over-intensified most of the highs and 
the lows. Since the effect of the non-adiabatic term on 
the average is one order of magnitude less than that of 
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Coef f . 


Case 


1 


Case 


2 


Case 


o 


heat- /cool- 
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| 
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+ 190 


1 
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j 

-33 
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4 193 


-16 
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* 

-23 
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1 

- 3 

J 


+ 197 


- 8 

i 


4 192 


25/20 


-13 


+ 270 


i + 4 


4 201 


1 0 


4 204 
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Table 3. Results of numerical verification for prognoses 
using advection and non-adiabatic terms with 
coefficients of 0.4 for warm advection and 0.5 
for cold advection. The coefficients used for 
the non-adiabatic term are shown in the table. 
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(Heating) 10 
^ (Cooling) 6 



Table 4. Results of numerical verification for prog- 
noses using all three terms of the tliichness- 
tendency equation, with coefficients as shown 
in the table. 
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the advection term, the advantage of including the non- 
adiabatic term would hardly appear. Secondly, for the non- 
adiabatic heat exchange, the zones of thermal concentration, 
where transformation of polar-continental air to polar- 
maritime air takes place, are really of primary concern. 
During the month of April, these zones of thermal concentra- 
tion were found along the northeast coasts of Asia and North 
America. For the three cases investigated, cyclonic activ- 
ity prevailed in these areas. With the 1000-mb lows al- 
ready over-intensified by the advection term and with the 
flow oriented such that heating took place in these areas, 
the lows were intensified further. For example, the low 
off the east coast of Japan in case 3, where the lowest 
observed height was -|— 30 feet, the prognosis with the ad- 
vection term alone gave a lowest height of - 170 feet while 
the non-adiabatic term decreased the lowest value further 
to - 250 feet. Thus, in spite of the fact that the behav- 
ior was in accord with theoretical expectations, the RMSE 
did not improve when the non-adiabatic term was included 
in this case. 

In order to investigate the effect of the non-adiabatic 
term further, a differential analysis was made by sub- 
tracting a prognostic map with both advection and non- 
adiabatic terms included (fig. 8) from the one with 
only the advection term. This isolated the effect of the 
non-adiabatic term (fig. 6) and it was found that by using 
a coefficient of 10 for heating, the largest 24-hr 

2 ? 



contributions by this term in areas of thermal concentra- 
tion were -282 feet off the east coast of Nova Scotia, at 
41N, 55W; and -230 feet off the west coast of Japan, at 40N, 
134W. 

Two possible sources of error warrant discussion. In 
the first place, the 1000-mb map was a relatively flat map, 
with the range of thickness values and 500-mb height much 
greater than the range of the 1000-mb height. Therefore, 
errors in the 500-mb height, when transmitted downward 
through the thickness, can cause serious deviations at the 
lower surface. Furthermore, it was obvious that one could 
not expect the locations of the pressure systems on the 
lower surface to be better than what the 500-mb prognosis 
can provide for this model. Secondly, the fact that the 
positions of the pressure systems were all better in the 
Atlantic Ocean than those in the Pacific Ocean, true in all 
three cases, was most probably due to the fact that weather 
reports in the Atlantic Ocean were more numerous and accu- 
rate than those in the Pacific Ocean. In other words the 
initial data from which the prognosis was started could 
have been in error to a certain degree in the Pacific Ocean. 
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5. Conclusions: 



Starting with a dynamical basis and utilizing the 
rapid computing capabilities of the CDC-1604 computer, 
empirical coefficients were found for each of the three 
terms of the prognostic equation designed to predict the 
1000/500-mb thickness. By subtracting the predicted 
thickness from the 500-mb 24-hr prognosticated stream 
function, a 24-hr prognosis of the 1000-mb map was obtained. 
The fact that this model did provide a map with most of the 
pressure systems appearing at the right places is encourag- 
ing. Due to the relatively small sample tested, the empir- 
ical coefficients found thus far were not necessarily the 
best ones; but they will nevertheless provide a starting 
point for further investigations of a similar nature. 

Improvements of the model might include the incorpor- 
ation of certain empirical rules in order to remedy the 
phenomena of over-predicting cyclogensis and anticyclo- 
genesis and the use of the weekly-mean sea-temperature 
field instead of the monthly mean. Furthermore, any im- 
provement in the prediction of the 500-mb maps will give 
increased accuracy at 1000 mb. In conclusion, it is the 
opinion of the present investigators that this model can 
provide a first approximation in the preparation of the 
surface prognostic map. 
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Figure 4. The observed map for 0000Z April 4 , 1955 . 
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Figure 5. The 24-hr prognostic map with all three 
terms of the prognostic equation in- 
cluded . 
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Figure 6. 



The 24-hr prognostic map with the thickness 
advection term only. 
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Figure 7. The 24-hr prognostic map with the 

thickness advection and the terrain 
terms. 
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•Figure 8. Hie 24-hr prognostic map with the thick- 
ness advection and the non-adiabatic 
terms. 
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Pigure 9. 



Analysis of the effects of the tefrair 
term. 
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Figure 10. 



The monthly-mean sea-temperature field 
of April. 
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